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ABSTRACT 

We use the Aquarius simulation series to study the imprint of assembly history on 
the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work 
regarding the influence of mergers on the mass density profile and the inside-out 
growth of halos. The inner regions that contain the visible galaxies are stable since 
early times and are significantly affected only by major mergers. Particles accreted 
diffusely or in minor mergers are found predominantly in the outskirts of halos. Our 
analysis reveals trends that run counter to current perceptions of hierarchical halo 
assembly. For example, major mergers (i.e. those with progenitor mass ratios greater 
than 1:10) contribute little to the total mass growth of a halo, on average less than 
20% for our six Aquarius halos. The bulk is contributed roughly equally by minor 
mergers and by "diffuse" material which is not resolved into individual objects. This 
is consistent with modeling based on excursion-set theory which suggests that about 
half of this diffuse material should not be part of a halo of any scale. Interestingly, 
the simulations themselves suggest that a significant fraction is not truly diffuse, since 
it was ejected from earlier halos by mergers prior to their joining the main system. 
The Aquarius simulations resolve halos to much lower mass scales than are expected 
to retain gas or form stars. These results thus confirm that most of the baryons from 
which visible galaxies form are accreted diffusely, rather than through mergers, and 
they suggest that only relatively rare major mergers will affect galaxy structure at 
later times. 

Key words: cosmology: dark matter - methods: N-body simulations - Galaxy : 
formation 



1 INTRODUCTION 

Hierarchical growth is a signature prediction of the ACDM 
cosmogony, our current standard picture of cosmic struc- 
ture formation. ACDM postulates a flat universe with a 
cosmological constant, cold dark matter, and gaussian ini- 
tial conditions generated at very early times. The basic 
units of nonlinear structure are dark matter halos that 
grow by accretion and merging as gas cools and condenses 
into galaxies in their cores. The statistics of this pro- 
cess are amenable to analytic modeling, which can in turn 
be validated and extended through cosmological N-body 
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The mass function of collapsed structure and its evolu- 
tion with time, the clustering of halos of different mass, the 
origin and character of scaling laws relating halo properties 
are all results that can be understood within the context 
of the excursion-set modeling making reference only to the 
initial power spectrum of density fluctuations and to the 
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universal expansion history (see, e.g., lLacev fc Cold 1 19931 ; 



iMo fc Whitell 19961 ; iNavarro et a"illl996h . 

The ACDM power spectrum, P(k), can be computed in 
detail using linear theory. Under the simplifying assumption 
that the "temperature" of the dark matter is zero (or, equiv- 
alenty, that the dark matter particle, if a thermal relic, has 
"infinite" mass), P(k) approaches k~ 3 on the smallest scales 
with a mass variance that diverges logarithmically there. 
Excursion-set theory then predicts that, at times of inter- 
est, effectively all of the mass of the Universe is in clumps of 
some mass. The assembly of a halo thus consists, in this sim- 
plified case, merely of the merging of the myriads of smaller 
mass subhalos that collapsed at earlier times. 

Under these conditions merging is the basic engine of 
halo growth. The rapid mixing driven by the violently fluc- 
tuating potential of a merger has awesome transformative 
powers. Mergers can erase, at least partially, memory of the 
initial conditions and leave remnants whose broad structure 
is roughl y independent of the cosmologic al conditions of for- 
mation jwhitei [197S1 ; Ivan Albadal llQSd '). Not all mergers, 
however, are created equal, and it has long been appreci- 
ated that the effects of major mergers differ qualitatively 
from those of minor events. In particular, major and mi- 
nor mergers affect differently the internal structure of the 
remnant. Major mergers lead to "rapid growth" in the mass 
of an object, which has been linked with radical changes 
in the halo structural parameters. Minor mergers, on the 
other hand, are associated with "slow growth" evolutionary 
phases that leave the inner structure of the main halo rela- 
tively intact and affect mainly the periphery of t he remnant 
ilSalyador-Sole et al.lll99Sl; IWechsler et al.ll2002l; IZhao et al.1 



120031 ; iTasitsiomi et al.ll2004l ; iDiemand et al.ll2007f ) 

The idea of merging as the exclusive mechanism of halo 
growth has received some backin g in the literature, most re- 
cently from iMadau et al. I (|200Sl ). who report that most of 
the mass in their Via Lactea simulation of a galaxy-sized 
halo "is acquired in resolved discrete clumps, with no evi- 
dence for significant smooth infall". On the other hand, a 
number of recent papers have also argued that the fraction 
of mass accreted "diff u sely" might be substan ti al (see, e.g., 
Fakhouri fc Mai l20ld ; lAngulo fc White! |2010| ; iGenel et all 



20101, and references therein). 



lAngulo fc White] (|2010l ). in particular, note that when 
a realistic cold dark matter particle candidate is chosen its 
small but non-negligible thermal velocity introduces a cutoff 
(and finite variance) in the power spectrum on small scales 
that can have a profound impact on the way the evolving 
hierarchy of collapsed structures develops. Working through 
the numbers appropriate for a neutralino-dominated Uni- 
verse, these authors argue that, as late as z ~ 20, most of 
the mass of the Universe is not yet part of any halo. These 
authors also argue that a typical galaxy-sized halo accretes 
at least 10% of its mass in diffuse form. 

If these numbers are correct the actual fraction of 
smoothly-accreted material in a typical N-body halo must 
be much higher, since simulations can only resolve a limited 
range of nonlinear scales and a fair fraction of the mass is 
expected to be locked u p in unresolved small mass clumps. 
lAngulo fc White! (|2010l ) argue that, even in the best sim- 
ulations currently available, up to 30-40% of the mass of 
a galactic halo could have been accr eted in diffuse form , in 
clear disagreement with the results of lMadau et~aH (|2008l ). It 




Figure 1. The evolution of the mass of the main friends-of- friends 
(FOF) progenitor of the six level-2 Aquarius halos. The curves 
show the mass in units of the mass at z = (labels on the left 
y-axis). Vertical segments indicate the mass ratio of the largest 
merger event occurring at each snapshot (labels on the right y- 
axis). Only merger events with mass ratio exceeding 0.05 are 
shown. Colors identify individual halos, as labelled in the figure. 
Note that only halo Aq-F-2 has undergone a major (/ m > 0.1) 
merger after z = 1. 



is clearly important to resolve this disagreement, especially 
given the importance of diffuse mass accretion fo r galaxy for- 
mation emphasized in recent papers (see, e.g., iKeres et al.l 
120051 ; iDekel et al.ll2009l . and references therein). 

We address these issues here u sing the N-body sim - 
ulations of the Aquarius Project (|Springel et al.l 20081 ) . 
This simulation series follows the formation of six different 
ACDM halos at various resolutions, and includes the best- 
resolved galactic dark matter halo simulated so far, an object 
with more than one billion particles within the virial radius. 
We begin with a brief description of the Aquarius Project 
in Sec.[2j and move on in Sec.|3]to a systematic study of the 
radial structure of halos in terms of the mass of their pro- 
genitor halos and the time of their accretion/merging. Sec. [3] 
considers the mode of accretion into the halos in detail and 
addresses the fraction of mass accreted in diffuse form. We 
end with a brief summary of our main conclusions in Sec.[S] 



2 THE NUMERICAL SIMULATIONS 

The Aquarius Project (|Springel et al.l l2008h consists of a 
suite of large N-body simulations of six dark matter ha- 
los of mass consistent with that expected for the halo of 
the Milky Way. Our simulations assume the ACDM cosmol- 
ogy, with parameter s consistent with the WMAP 1-year data 
(|Spergel et al.ll2003h : matter density parameter, Qm = 0.25; 
cosmological constant term, Qa = 0.75; power spectrum nor- 
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Figure 2. Mass accretion history of halo Aq-A. Curves in black 
correspond to the level-2 resolution halo, in red to the level-4 run. 
All masses are normalized to the FOF mass of the halo at z = 0. 
Four definitions of halo mass buildup are compared. The long- 
dashed (top) curves indicate the mass of all particles "associated" 
with the halo; i.e., those that were part of the main progenitor at 
any time before redshift z. The dotted (bottom) curves are the 
mass of particles that belong to the most massive progenitor at 
all times after z (i.e., after the time of last accretion). The thin 
solid curves (second set from top) indicate the cumulative mass 
of particles as a function of the redshift of first accretion into the 
main progenitor (regardless of whether they leave and re-enter the 
main progenitor subsequently). The short-dashed curves indicate 
the conventional FOF mass of the main progenitor at each time. 
This comparison illustrates the fact that a significant fraction 
of the mass of the main progenitor is pushed out of the halo 
boundary during its evolution; much of it is re-accreted later, but 
more than 20% is still outside the main halo at z = 0. 



malisation, erg = 0.9; spectral slope, n s — 1; and Hubble 
parameter, h — 0.73. 

The halos were identified in a 900' ! -particle N-body 
simulation of a cubic volume 100 h -1 Mpc on a side, a 
lower resolution version of the Millennium-II Simulation 
jBovlan-Kolchin et aI.ll2009T ). This volume was resimulated 
using exactly the same power spectrum and phases of 
the original simulation, but with additional high-frequency 
waves added to regions surrounding the initial Lagrangian 
volume of each halo. The high-resolution region was popu- 
lated with low- mass particles and the rest of the volume with 
particles of higher mass. These "zoomed-in" simulations of 
selected regions or individual objects have become common 
practice to make g alaxies; for details we refer the reader to 
iPower et al.1 (|2005t ). 

The six Aquarius halos are labelled "Aq-A" through 
"Aq-F" . Each was resimulated at different resolutions in or- 
der to assess numerical convergence. A suffix, 1 to 5, iden- 
tifies the resolution level, with level 1 denoting the highest 



resolution. Between levels 1 and 5, the particle mass ranges 
from m p ~2x 10 3 Mq to ~ 3 x 1O 6 M . Level 1 was per- 
formed only for Aq-A and contains roughly 1.1 billion parti- 
cles within the virial radius. All six halos were simulated at 
level-2 resolution. Each of these has more than 100 million 
particles within the virial radius. 

In this study, we analyze primarily the level-2 simula- 
tions but we also use lower resolution versions of Aq-A to 
test for numerical convergence. For the level-2 simulations, 
the particle mass is m p ~ 1 x 10 4 h^ 1 Mq and the softening 
length is e = 48ft _1 pc. At z = 0, the six haloes have a sim- 
ilar "virial" mass, M 2 oo ~ 1-2 x 10 12 h^ 1 Mq, where M 20 o 
is the mass contained within r2oo, the radius of a sphere 
of mean density 200 times the critical density for closur^E 
The circular velocity curve of the halos peaks at roughly 
V m ax = 220 ± 40 km s _1 . For fur ther details of t he Aq uar- 
ius Project, we refer the reader to lSpringel et al] (|2008l ) and 
iNavarro et ail (|2010T ). 

At every snapshot in the simulation we find nonlinear 
structures using t he friends-of-friends (FOF) algorithm of 
iDavis et al.l (|l985l ). with a linking length of 0.2 times the 
mean interparticle separation and 32 particles as the min- 
imum number of particles per group. We then construct a 
merger tree for the final Aquarius halos linking FOF pro- 
genitors at each time. We also identify bound substructures 
within e ach FOF halo ( subha los) using the SUBFIND algo- 
rithm of lSpringel et al.l d200of). Merger tree s for su bhalos are 
constructed as described in lspringel et all (|2008l ). 

For simplicity, unless otherwise explicitly noted, we 
shall identify a halo with the FOF structure that contains it. 
Note that the mass of FOF halos does not necessarily coin- 
cide with the virial mass alluded to above (FOF structures 
are larger; they typically enclose a halo and a small part 
of its surroundings); we shall comment on these differences 
when appropriate in the analysis that follows. 

Fig. Q] shows the growth of the FOF mass of the main 
progenitor of each halo, normalized to its value at the 
present time, z — 0. Merger events with mass ratio greater 
than fm = Mp r0 g/M ma i n = 0.05 are noted by vertical lines 
in the colour corresponding to the appropriate halo (scale 
on the right). 

If we define the formation time of a halo as the time 
when the main progenitor first reaches half the final mass, 
then the formation time of the halos is z ~ 1.2 — 2.2, except 
for halo Aq-F-2 which is clearly different from the rest. Its 
formation redshift is z = 0.6, when its mass almost dou- 
bles as a result of an almost equal-mass merger (the actual 
mass ratio of the two progenitors is f m = 0.75). The other 
five halos have similar mass growth histories but different 
merger histories. For example, halos B, D, and E experi- 
enced major mergers at high redshift (z ~ 2 — 7), while 
halos A and C grew in a relatively quiescent fashion and did 
not experience any major mergers after z — 6. Although the 
Aquarius haloes have similar final masses, they ha v e var - 
ied formation histories, which iBovlan-Kolchin et al.1 (|2010T l 
has shown, sample the range of behaviours seen in the Mil- 
lennium II simulation for halos of this mass. The Aquarius 
haloes therefore provide a suitable sample to study the diver- 



1 This choice defines implicitly the virial radius of the halo, r200; 
and its virial velocity, V200- 
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0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 

r [kpc/h] M pfog r [kpc/h] M pfog r [kpc/h] M pfog 

Figure 3. The radial distribution of particles in the FOF 2 = halo colour-coded according to M pl - og , the mass of the progenitor to 
which each particle belonged at the time of (first) accretion into the main halo. The bars represent the fraction of the mass in each 
spherical shell brought in by halos with mass in the range indicated by the key to the right of each panel. This key also gives the total 
fraction (summed over all radial shells) of mass brought in by different progenitors. Material accreted before z = 6 is indicated in grey. 
Diffuse material, that is, particles that were not part of any FOF halo at the time of accretion, are indicated in red. Masses are in units 
of h,- 1 M . The white curve gives the fraction of the total FOF halo mass in each radial shell. 



sity in assembly histories of halos similar to that surrounding 
the Milky Way. 



3 MASS AND ACCRETION HISTORY 

3.1 Mass growth and definition of accretion 

Starting at z = we identify the main trunk of the merger 
tree of each of our halos by stepping back in time, defining 
the main progenitor at time n — 1 to be the largest FOF 
halo which is a progenitor of the main progenitor at time n. 
For each particle in the final object at z = 0, we register the 
redshift when it was accreted into the main progenitor (i.e. 
when it first ceased to be a single particle or a member of 



( Balogh et al 



another FOF group), 



and the mass of the FOF group, 



Mp r0 g, to which it belonged at the snapshot immediately 
preceding the time of accretion. For "diffuse" or "smooth" 
accretion, terms we use interchangeably throughout, M pros 
equals the mass of a single particle. 

This seemingly straightforward definition of accretion 
is complicated by the fact that some particles can leave the 
main progenitor and be re-accreted again later on. This pro- 
cess can actually recur multiple times, and is usually associ- 
ated with accretion events, where a small but non-negligible 
fraction of the mass is propelled into highly energetic orbits 



iLudlow et al.l 



2000; I Gill et all 120051 : iDiemand et al ] 120071 : 



20091 ). 



There is therefore some ambiguity in the meaning of 
accretion time whose effects we illustrate in Fig. [2] Here we 
compare the growth of halo Aq-A using several plausible 
definitions of accretion and two different levels of resolution 
to check for possible numerical artifacts. The dashed curves 
(second from the bottom) track, as in Fig. [T] the conven- 
tional FOF mass of the main progenitor, normalized to its 
value at z = 0. The long-dashed curves (top) show, on the 
other hand, the cumulative mass of all particles that, at any 
time before 2, have been part of the main progenitor. These 
"associated" particles exceed the FOF halo at z = by 
more than 20%, highlighting the importance of the energy 
redistribution process described in the previous paragraph. 

The solid thin curves, on the other hand, track the 2 = 
FOF particles, but use the time of first accretion to define 2. 
The difference between this and the conventional FOF mass 
is a direct indicator of the accretion-escape-reaccretion pro- 
cess alluded to above. Finally, the dotted (bottom) curves 
use the time of last accretion of particles in the 2 = FOF 
group. Fig. [2] shows clearly that a halo is a dynamic object 
and not a static "bucket" of mass that gets progressively 
filled by accretion. For simplicity we shall in what follows 
adopt the time of first entry as our default definition of ac- 
cretion but we caution that other definitions of accretion 
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Figure 4. As Fig [3] but for z acc , the redshift of (first) accretion of particles into the main halo. Different colours indicate accretion 
redshift in the intervals indicated by the key to the right of each panel. This key also gives the total mass fraction (summed over all 
radial shells) accreted in each redshift interval. The white curve gives the fraction of the total FOF halo mass in each radial shell. 



may on occasion be more useful, depending on the aim of 
the analysis. The good agreement between the results for 
Aq-A-4 and Aq-A-2 show that these conclusions are insen- 
sitive to the numerical resolution of the simulations. 



3.2 Progenitor mass distribution as a function of 
radius 

We now investigate the fate of particles accreted in pro- 
genitor halos of different mass. More specifically, we investi- 
gate the distribution of halo mass in spherical shells (centred 
at the location of the minimum in the gravitational poten- 
tial) and apportion the contribution according to the mass, 
Afprog, of the subhalos that brought each particle into the 
main progenitor. 

This is shown in Fig. [3] which gives, in each panel, the 
fraction of particles in concentric radial shells of each z = 
Aquarius halo, split into six mass bins, according to the mass 
of the progenitor halo at the time of accretion. Each bin is 
identified by a different colour according to the key to the 
right of each panel. This key, in addition, gives the total 
fraction (summed over all radial shells) of mass brought in 
by progenitors of different mass. Material that was already 
in place at z — 6 (when few resolved progenitor halos exist) 
is indicated in grey; in this section and the next we do not 
consider the assembly history at earlier times. Particles that 
in the snapshot preceding accretion are unattached to any 
resolved halo (i.e., diffuse accretion) are indicated in red. 



Fig.0sb.ows that there is considerable halo-to-halo vari- 
ation in the mass spectrum of the progenitors of the final 
halo. Consider, for example, the integrated halo mass. As the 
key on the right of each panel shows, the fraction of mass 
that was already in place at z = 6 ranges from ~ 5 to ~ 15%. 
The fraction of mass accreted diffusely is substantial in all 
cases, ranging from ~ 30 to ~ 40%, The diffuse fraction in- 
creases with radius, from a few per cent within 10 h~ x kpc, 
to more than 20 per cent within 100 h~ x kpc . This agrees 
with the results found by iHelmi et alj (|2002l ) from an N- 
body simulation of a cluster halo scaled to a galactic mass. 
The radial behaviour is diverse. For example, in halos A and 
C the material in the innermost region, r < 1 kpc, was al- 
ready in place before z = 6 and has undergone little change 
since. By contrast, in halos B, E and F, most of the central 
mass was brought in after z — 6 through mergers involving 
host halos with mass greater than 10 h~ Mq. 



Halo D is intermediate between these two extremes. 
Large progenitors, of mass > 10 10 h -1 Mq, bring in between 
20 and 40% of the final mass at all radii. In general, the 
larger the mass of the carrier halo, the greater the proba- 
bility of the particles ending up in the central regions. Very 
little of the diffuse material makes it into the central regions, 
r < 10 ft -1 kpc; indeed, most of it stays in the outer regions 
where it typically contributes most of the mass. 
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Figure 5. The mean accretion time of particles in different radial 
shells at z = 0. Accretion times are shown relative to the forma- 
tion time of each halo, if orm , defined as the time when the halo 
first reached half its final mass. When defined this way the mean 
accretion time profile is similar for all Aquarius halos. 



3.3 Accretion time distribution as a function of 
radius 

We now turn our attention to the radial distribution of par- 
ticles as a function of their accretion time. Analogously to 
Fig. OH Fig. |4] shows the distribution of accretion redshifts, 
binned according to distance to the centre of the halo. In 
each radial shell, the contribution from material accreted in 
different time intervals is shown by strips coloured according 
to the key shown to the right of each panel. The key also 
gives the fractions of the total halo mass (summed over all 
shells) accreted in each redshift interval. 

The inside-out nature of halo assembly is clearly appar- 
ent in Fig. [4] On average, the peak contribution from each 
accretion redshift interval marches outwards with time. The 
inner regions are populated mostly by particles that were 
accreted early; the outer layers were added gradually later. 
The cores of halo A and C were in place before z = 6 and 
evolved little thereafter. In halos D, E, and F, the core par- 
ticles were accreted by z = 3, but for halo B, the core is 
accreted at z = 1 because two major mergers in the redshift 
interval 1 < z < 3 bring in almost 50% of the core mass (see 
Fig- [!])■ These mergers happen relatively early, while the to- 
tal halo mass is small, and disrupt the original core which 
then reforms from the new material. By contrast, the late 
major merger undergone by halo F has a relatively minor 
effect on the core, probably because of the orbital param- 
eters of the merger. The core in this case is actually made 
primarily out of material that was accreted in earlier, lesser 
mergers at z ~ 4. 

The radial dependence of accretion time is quantified 




0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 

r [kpc/h] M pro , r [kpc/h] M prog 

Figure 6. A convergence test of the radial gradients in progenitor 
mass and accretion redshift. The top two panels are analogous to 
Fig. [4] the bottom panels to Fig. \3\ but for the level-2 and level- 
4 Aq-A runs. These runs differ solely in numerical resolution; 
level-2 runs have 30 X more particles and 5 X smaller gravitational 
softening than their level-4 counterparts. The excellent agreement 
shows that the results presented in Figs. [3] and [4] are insensitive 
to numerical resolution. 



further in Fig. [5] which shows the average accretion time as 
a function of radius. Accretion time is plotted relative to the 
formation time of the halo, iform, defined as the time when 
the main progenitor first reaches half its final mass. With 
this normalization, the radial dependence of the accretion 
time is fairly similar for all Aquarius halos. On average, the 
material in the inner 10 h~ kpc is assembled 2 to 4 Gyr 
before tf orm whereas the material beyond 100 ft -1 kpc falls 
in 2 to 4 Gyr after tform- This onion- like growth is generic for 
cold dark matter halos of galactic sca le; it was seen also in 
a scaled cluster N-body simulation bv lHelmi et al.l l|2003l ). 

3.4 Numerical convergence 

Before discussing these results further we should verify that 
the trends presented above are not unduly influenced by 
numerical resolution. The availability of simulations of the 
same halo at varying resolution allows for direct testing of 
the reliability of our results. We do this by comparing the 
level-2 simulation of halo Aq-A, which is the one we have 
analyzed so far, with its level-4 counterpart. The level-4 sim- 
ulation has about 30 x poorer mass resolution and 5x poorer 
spatial resolution (softening). 

The test is carried out in Fig. [6] which shows the distri- 
butions of accretion time and the mass spectrum of progeni- 
tor halos as a function of radius. These figures are analogous 
to Figs. [3] and [4] Panels on the left show the level-4 results, 
those on the right the level-2 results. It is clear that the 
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Figure 7. The fraction of particles in a series of spherical shells 
that were accreted smoothly (black circles - solid line), by minor 
mergers (starred symbols - dotted) and by major mergers (open 
triangles - dashed). The distinction between minor and major 
mergers is made at a mass ratio of 10:1. Fractions are averaged 
over all six level-2 Aquarius halos and error bars show the rms 
scatter among halos. Despite substantial scatter, the trends are 
clear. Major mergers contribute to the inner regions, diffuse accre- 
tion and minor mergers mainly to the outer regions which contain 
the bulk of the mass. 



convergence of these properties is excellent. There is no dis- 
cernible difference in the distributions of z a cc and at most a 
~ 10% difference in the distributions of M prog in the lowest 
mass range, M prog < 10 7 h^ 1 M©. 



4.1 Major mergers vs minor mergers 

In Fig. [7] we illustrate how major mergers, minor mergers 
and diffuse accretion contribute to the z = mass in a series 
of spherical shells, each spanning a factor of two in radius. 
The symbols joined by lines give results averaged across the 
six Aquarius halos, while the error bars indicate the rms 
scatter among halos. Within ~ 10 h^ 1 kpc, major mergers 
are the dominant source of the material, providing typically 
40% of the mass, while minor mergers and smooth accretion 
bring in about 30% each respectively; within ~ 1.0 h~ x kpc, 
major mergers contribute more than two thirds of the mass. 
Note, however, that less than 10% of halo mass lies within 
~ 10 h^ 1 kpc and less than 1% lies within ~ 1.0 /i -1 kpc. 
Note also from Fig.|4]that the great majority of these major 
mergers occurred at z > 3 and many of them at z > 6. Only 
in halos B and F are there substantial contributions to these 
regions from major mergers at redshifts below 3. The large 
error bars on these points indicate that the scatter of the 
major merger contribution to the inner regions of halos is 
large. 

Beyond 10 /i" 1 kpc, in the region which contains the 
bulk of the halo mass, both minor mergers and diffuse ac- 
cretion contribute more to halo growth than major mergers. 
Indeed, averaged over all six halos, major mergers contribute 
only 17% of the total mass growth, with the values for indi- 
vidual halos ranging from 3% (Aq-A) to 36% (Aq-F). For a 
stricter definition of a major merger, requiring a mass ratio 
of 3:1 or less, the mean major merger contribution drops to 
just 9%, with individual values ranging from < 0.1% (Aq- 
A, Aq-C) to 25% (Aq-F). Thus, major mergers are typi- 
cally a small contribution to overall halo growth. The rest 
is split almost evenly between minor mergers and "diffuse" 
accretion. It is interesting that the scatter in each of these 
contributions is very close to half that in the major merger 
contribution. This shows that the minor merger and diffuse 
fractions fluctuate up and down together, with minor merg- 
ers contributing slightly less than half of the material not 
accounted for by major mergers at each radius and in each 
halo. 



4 MODES OF ACCRETION 

In this section we study how the growth of halos is appor- 
tioned between major mergers, minor mergers, and smooth 
accretion; how the material added in these modes is dis- 
tributed in radius in the final halos; and how much vari- 
ation there is between halos. We will adopt a FOF mass 
ratio of 10:1 as our standard division between major and 
minor mergers, although we will also give some results for 
the stricter 3:1 ratio adopted as a boundary by some au- 
thors. The fact that we limit our FOF group catalogues to 
systems with at least 32 particles means that the boundary 
between minor mergers and "smooth" accretion occurs at 
a mass ratio of about 10 :1 at z = dropping to about 
10 5 ' 5 :1 at z = 4 and to even smaller values at higher red- 
shifts. In this section we consider increases in mass through 
each of these growth modes throughout the entire history 
of each halo, rather than halting at z = 6 as in previous 
sections. 



4.2 Diffuse accretion 

Given the conflicting claims in the literature regarding the 
importance of diffuse accretion discussed in the Introduc- 
tion, it is important to explore possible biases and sub- 
tleties involved in reckoning the amount of diffuse mass 
accreted. The dynamic nature of halo buildup highlighted 
above (Sec. 13. 1[) introduces amibiguities in the meaning 
of accretion, so we compare four alternative definitions of 
/smooth, the total fraction of mass in the FOF halo at z = 
that has been added smoothly: 

I: all particles that were not part of any 32+ particle 
FOF group in the snapshot immediately before the time of 
first accretion, z acc ; 

II: all particles that were never part of any 32+ particle 
bound structure (as identified by SUBFIND) before z acc ; 

III. same as (II) but for 20+ particles; 

IV: same as (I) but for all snapshots before z a cc- 

Fig. [S] compares results for the six level-2 Aquarius ha- 
los. Criterion (I), probably the simplest, is seen to give the 
largest estimate of / sm0 oth in all cases. This criterion omits 
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those particles that were part of FOF halos in the past, but 
that have left them and are unattached to any resolved struc- 
ture just before accretion. These make a surprisingly large 
fraction (about half\) of the smooth fraction computed using 
criterion (I), as shown by the bottom curve corresponding 
to criterion (IV). 

One shortcoming of criterion (IV), however, is the pos- 
sibility that FOF groups may artificially link in physically 
unrelated particles. This is especi ally true in small -N groups 
(for a recent discussion, see, e.g.. iBett e t al. 200?]). Criteria 
(II) and (III) account for this by requiring particles to be 
part of bound structures; varying the threshold from 20 to 
32 particles has negligible effect on the results. This extra 
condition is seen to increase / am0 oth by roughly 50% rela- 
tive to criterion (IV). The contribution of diffuse accretion 
seems, therefore, to be genuinely high, between 20 and 40% 
of the final halo mass overall. 

We note that, strictly speaking, /smooth depends on the 
total number of snapshots used in its estimation. The num- 
bers quoted above are based on a total of 128 snapshots, 
but for one of the runs (Aq-A-2) data were stored for 1024 
snapshots. The estimate of /smooth according to criterion (I) 
changes little when considering 1024 or 128 snapshots: from 
41% to just 36%. The changes are even smaller for defini- 
tions II or III. 

Finally, we consider the dependence of /smooth on the 
mass resolution of the simulations. For this, we use four dif- 
ferent realizations of halo Aq-A, from level 2 to level 5. The 
results are shown in Fig. [9] (criteria II and IV: open diamonds 
and asterisks, respectively). As expected, there is a system- 
atic decrease in /smooth with increasing resolution, measured 
in Fig. [5] by M cu t, the mass of a group of 20 particles. 

Given that / sm ooth depends on resolution, we need to 
ask how secure is our estimate of this quantity from the 
simulations. We can answer this by analyzing Monte Carlo 
merger trees built using the excursion-set formalism, con- 
strained at z = to make a halo of mass comparable to those 
in the Aquarius s et. In parti cular, we u s e the algorithms of 
iParkinson et al. | (|2008h and ICole et al. I (|2008h . which were 
tuned to m atch the N-body mer ger trees of the Millennium 
Simulation ()Springel et al.|[2005l l . This approach has the ad- 
vantage that a cutoff mass can be easily introduced in order 
to mimic the limited resolution of a simulation (see, e.g., 
lAngulo fc Whitell2010D . 

The excursion-set results are shown by the dotted line 
in Fig. [8] and by the connected dots in Fig. [9] The Monte 
Carlo trees, when trimmed to match the resolution of the N- 
body simulations, give results in good agreement with the 
simulations. The theoretical calculation also confirms the 
large scatter in /smooth seen in the simulations. (The "error 
bars" on the Monte Carlo results denote the rms scatter 
among several hundred realizations.) 

The trends shown in Fig. pimply that further improve- 
ments in resolution would result in only small reductions in 
the value of /smooth- Indeed, /smooth seems to depend more 
strongly on the particular definition adopted for smooth ac- 
cretion than on numerical resolution, at least for the 100- 
million particle halos we consider here. The vertical line in 
Fig. [9] shows the value of M cut corresponding to Aq-A-1, 
the best-resolved, billion-particle halo in the Aquarius set. 
The Monte Carlo tree results suggest that its additional res- 




C D E 
Haloes 



Figure 8. The fraction of mass accreted in diffuse form in each 
level-2 Aquarius halo. The x-axis lists the name of each halo. 
Each curve corresponds to one of the definitions of "diffuse ac- 
cretion" introduced in the text. Briefly, (I) are particles that are 
unattached to any FOF group identified in the snapshot imme- 
diately before first accretion; (II) refers to particles that do not 
belong to any bound structure with ^ 32 members in any snap- 
shot before first accretion; (III) is as (II) but for 20 members; and 
(IV) denotes material that did not belong to any N ^ 32 FOF 
group at any time before first accretion. All curves use a total of 
128 snapshots to estimate /smooth- The excursion-set prediction 
for a halo of the same mass and comparable numerical resolution 
is shown by the horizontal dotted line. 



olution would result in only a very small decrease in the 
smoothly accreted fractionQ 

We conclude from this exercise that the substantial frac- 
tion of mass found to be accreted smoothly in our simula- 
tions is a robust result. We can ther efore confidently rule 
out the claim by iMadau et alj (|2008h that at most 3% of 
the mass of a galaxy-sized halo can be supplied by smooth 
accretion. It is not clear at this point what the cause of the 
disagreement is, but it is likely to be related to the way in 
which these authors compute diffuse accretion rather than 
to differences in the simulations themselves. For example, 
the mass they regard as having been accreted in "identifi- 
able subunits" is just the sum of the masses of all progenitor 
halos that contain particles that make it into the final sys- 
tem. A substantial fraction of that summed mass includes 
particles that are not part o f the halo at z — 0; this could 
have led lMadau et al.1 (|2008h to overestimate the mass con- 
tributed by discrete identifiable subunits and, therefore, to 
underestimate / sm0 oth- 



2 We have not attempted the analysis presented here for Aq-A-1 
because of the formidable computational task involved in building 



full merger trees for this simulation. 
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Figure 9. The dependence of the diffuse accretion fraction, 
/smooths on the mass resolution of the simulations, M cu t, defined 
as the mass of a group of 20 particles. Open diamonds and aster- 
isks refer to diffuse accretion definitions II and IV, as described in 
the text and in the caption to Fig. [8] Connected symbols refer to 
resolution levels 2 through 5 for halo Aq-A. Colors indicate differ- 
ent halos, as labelled in the panel. Dots with error bars indicate 
the mean and rms scatter in several hundred Monte Carlo assem- 
bly histories constructed from excursion-set theory. The vertical 
dotted line indicates the value of M cu t for Aq-A-1, the highest 
resolution simulation in the Aquarius series. 



Our conclusions are in a greement with t hose of 
lAngulo fc White] (|2010l l (see also lGenel et alj|2010h : merg- 
ers and smooth accretion are both defining features of the 
hierarchical buildup of a CDM halo. 

At redshifts z ^ 6 the intergalactic medium is fully 
photoionized and gas is unable to collect in halos with max- 
imum circular velocity below about 15 km/s, correspond - 
ing to masses below roug hly 1O 8 M jOkamoto et aI.ll2008T ). 
This is well above the resolution limit of our simulations, so 
the gas associated with these low-mass halos should be con- 
sidered to be diffusely accreted along with that associated 
with the "diffuse" dark matter. Taking the limit at exactly 
10 8 /i _1 M Q for simplicity^, we find that, on average, our ha- 
los accrete 56% of their baryons diffusely, with the numbers 
for individual halos ranging from 46% (Aq-F) to 64% (Aq- 
C). The bulk of baryonic accretion is thus predicted to be 
smooth rather than clumpy for objects of Milky Way scale. 



4.3 Evolution of the inner mass profile 

As Fig. [7] indicates, major mergers contribute, on average, 
just under half of the particles in the inner 10fe _1 kpc of 
Galactic halos. This is the region occupied by the luminous 



The following numbers are insensitive to this choice. 




0.0 



2 4 6 8 10 

Age of the Universe ( Gyr) 

Figure 10. Total mass enclosed within different physical radii, 
A/(< r), for r = 1; 2; 8; and 32 kpc, as a function of cosmic 
time. Different colours correspond to different halos, as labelled 
in the figure. Masses in each panel are normalized to their values 
at z = 0. Time is labelled at the bottom, redshift at the top. 
Except for halo Aq-F, which is the remnant of a recent major 
merger, the inner mass profile of Aquarius halos has been very 
stable for the past 5-6 Gyr, a period comparable to the age of the 
Solar System. 



component of the central galaxy, and it is thus interesting to 
analyze in detail how the mass profile in this region evolves 
with time. A thin stellar disk, for example, could react to 
clumpy addition of material by thickening and becoming 
dyamically hotter, potentially violating obser vations of the 
thin disk in the solar neighbourhood (see e.g. iBenson et al.l 
l2004h . 

We investigate the stability of the inner halo explicitly 
in Fig. 1101 where we plot the mass enclosed within 1, 2, 8, 
and 32 (physical) kpc0 from the centre of each halo as a 
function of time. Each curve is normalized to the enclosed 
mass at z = 0. Except for halo F, which undergoes a major 
merger at z ~ 0.6, all the halos show exceptionally stable 
inner mass profiles over at least the past 5 Gyr (z < 0.6), 
the age of the Solar System. Five out of six Aquarius halos 
could, in principle, host a disk as thin and cold as that of 
the Milky Way. Late-accreted mass typically settles in the 



4 We assume h 
the Milky Way. 



0.73 in order to compare with observations of 
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outskirts of a halo, thus allowing the hierarchical growth of 
halos to be reconciled with the ubiquitous presence of thin 
stellar disks. 



5 SUMMARY AND CONCLUSIONS 

We have analyzed the build-up of six ACDM halos simu- 
lated as part of the Aquarius Project to study the influence 
of assembly history on halo structure. We focus on the distri- 
butions of progenitor mass and accretion time for particles 
at different radii at the final time, and discuss various plau- 
sible definitions of accretion time, together with the difficul- 
ties involved in estimating the total mass fraction accreted 
smoothly. We compare simulations of the same halo carried 
out with different resolution in order to assess the sensitivity 
of our results to numerical limitations. 

Although there is considerable variation from halo to 
halo, our simulations exhibit a number of very clear trends. 
Our main conclusions may be summarized as follows. 

• There is a strong radial gradient in accretion time, 
which confirms that halos are built from the inside out. 
Later accreting material settles farther from the centre of 
the halo; particles that today reside inside 10 /i -1 kpc are 
typically accreted ~ 3 Gyr earlier than particles that reside 
at 100 ft -1 kpc from the centre. 

• Similarly strong correlations exist between distance of 
a particle from halo centre and the mass at accretion of the 
progenitor halo which contained it. The innermost regions 
are dominated by particles brought in by massive clumps, 
"major mergers" with mass ratios exceeding 1:10, as well 
as by those that joined the main progenitor at very early 
times (z > 6). Mass accreted diffusely and in minor mergers 
predominantly populates the more distant parts of the halo 
and dominates the total mass. 

• Minor mergers and diffuse accretion contribute approx- 
imately equally to the mass of each halo at each radius, at 
least at the resolution of our six Aquarius halos where mi- 
nor mergers can be distinguished up to mass ratios exceeding 
10 6 :1. 

• The inner mass profile of a halo is very stable at late 
times in systems that stay clear of major mergers. In five 
of the six Aquarius halos the profile within 32 kpc barely 
changes in the past 5-6 Gyr. 

• Diffuse accretion contributes a substantive fraction of 
the final mass of the halo, roughly 30-40% in our simulations. 
This is a robust result compatible with expectations from 
excursion-set modeling. 

• Our analysis shows that some of the material accreted 
smoothly had previously been part of other collapsed struc- 
tures, from which it was probably ejected by mergers. The 
same mechanism leads a fair fraction of particles in the main 
halo to cycle in and out of its main progenitor; at z = the 
cumulative mass of all particles "associated" in the past with 
the main progenitor exceeds the final mass of the halo by at 
least 20%. 

• After reionization, gas is unable to collect in dark halos 
with masses lower than about 10 8 ft -1 Mq. As a result more 
than half of all the baryons associated with halos of Milky 
Way scale are expected to be accreted smoothly, rather than 
in clumps. 



These results emphasize the dynamic nature of halo 
buildup and provide insight into the radial structure of a 
halo and the history of its assembly process. The view that 
emerges highlights some misconceptions regarding hierarchi- 
cal growth. CDM halos are not passive repositories where 
mass is added in continuous but discrete events, but rather 
lively systems that can lose as well as gain material through- 
out their lifetimes. Diffuse accretion, recurring infall, escape 
and fallback, are all processes that play important roles in 
the build-up of CDM halos. 
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